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INTRODUCTION 


This report presents a very fast hidden line technique (JONES-D) developed by 
the author, for the interactive graphical display of NASTRAH finite element 
models. It is the author's opinion that hidden line plotting together with 
haloed line and normal (all lines visible) plotting is required for an 
effective finite element plot package. Illustrated in figure 1 are hidden 
line, haloed line, and normal plots of the same finite element model. In this 
example, the hidden line plot is clearly the most effective form for human 
understanding. Haloed plotting is the most effective with models containing 
no, or very few, surfaces; thereby rendering hidden line plotting 
inappropriate. Techniques to perform hidden line plotting have been much 
discussed beginning with the advent of computer graphics in the eaily 1960* s 
and continuing into the present era. Given the bulk of this prior work (for 
example references 1 through 9), why develop a new method? The answer is that 
these prior methods appear to lack the speed for effective interactive use, or 
lack the features required to effectively plot NASTRAN models. Except for the 
Watkins technique (reference 9), code to implement these techniques were not 
published or generally available. Experience in using the Watkins technique 
had shown it not to be acceptable for the proposed use. References 10, 12, 
and 13 were published after the JONES-D method had been completed or 
substantially so, these new methods together with the Watkins method are 
compared with the technique developed herein in a subsequent section of this 
report. 


DLSIGN GOALS 


In order to understand the design goals, the computing environment must be 
defined. The Applied Engineering Division's (NASA Goddard Space Flight 
Center) computer system consists of a Digital Equipment Corporation VAX 11/780 
computer and support devices. The host operating system used for developing 
and testing the hidden line routine was VAX VMS version 2.5. The maximum 
interactive user working set size was 256 kbytes. The target graphic display 
devices were the four Tektronix 4014 (or equivalent) terminals used by the 
structural analysts within the Applied Engineering Division. These terminals 
were connected to the host VAX via 9600 baud RS-232 data links, and the base 
plot software was Tektronix PLOT 10 running on the host. The specific design 
goals for the hidden line routine were: 

a. The hidden line routine must function within the NPLOT NASTRAN plot 

package being developed by the author. The NPLOT program supports plotting of 
normally used NASTRAN elements: line elements (CBAR, CROD, etc.), surface 

elements (CQUAD, CTRIA, etc.) and solid elements (CHEXA and CPENTA). The 
hidden line technique must support these elements. The hidden line routine 
must also be supportive of labeling of visible grid points and elements. 

b. Hidden line plotting must be as responsiye to the interactive user as 
normal, all lines visible, plotting. This requires an elapsed time per hidden 
line plot of about the same as, or less than, that for a normal plot. It was 
recognized that this was a very ambitious goal. If we examine the data flows 
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associated with the normal, halo, and hidden line plotting, figure 2, we can 
see why this goal is at least conceptually possible and why halo line plotting 
is usually less responsive than normal plotting. The data flow for normal 
plotting consists of processing N vectors through PLOT 10 software and 
transmitting the appropriate conmands to the graphic terminal. The halo plot 
module, figure 2, normally expands the N vector list into a vector list of 
length l.b N. Therefore, we would expect halo plotting to be more than 1.5 
times slower than normal plotting. Conversely, the hidden line ntodule, figure 
2, truncates the input vector list of length N to about 0.5N. Thus, hidden 
line plotting would equal the responsiveness of normal plotting, provided that 
the hidden line module can execute in about the same time as that for 
processing 0.5 N vectors through the normal plot flow path. 

c. A high degree of reliability is required. A plot routine that 
typically generates plot errors is worse than useless for the debugging of 
finite elenient models. The required error rate from the hidden line routine 
must be very low. 

IMPLEMENTATION 


Several different variations of the same basic hidden line method were 
sequentially developed in the course of this effort. To keep track of the 
different versions, they were assigned the names JONES-A through JONES-D. The 
fastest and niost recent version, JONES-D, is the subject of this report. The 
algorithni (Appendix B) was written in DEC VAX FORTRAN 77 and embedded in the 


NPLOT NASTRAN plot program. This implementation made considerable use of the 
VAX virtual memor> feature to maximize performance. The basic flow of the 
technique is as follows: 

0 INPUT: The chief inputs to JONES-D from NPLOT are the global edge 
vector list, global surface list and grid point table. It should be noted 
that NPLOT operates to produce nonredundant global edge vector and surface 
lists. 

0 PREPARATION: The vector and surface lists are operated on to produce 

vector and surface data arrays to speed the down stream computations. For 
example, the minimum/maximum values for each vector and surface in the display 
co-ordinate system is computed. Spatial sorting of the vector and surface 
data is performed. Illustrated in figure 3 is a simplistic view of this 
technique. Based on the complexity of the model, n X n mesh X-Y grids are 
imposed on the model and lists of vectors and surfaces are generated for each 
spatial cell via bucket sorting. Separate grids are used for the vector and 
surface X-Y bucket sorts. The grid density for the vector X-Y sort is based 
on the number of edge vectors due to surface or solid elements; for surface 
X-Y sorting the grid density is based on the number of surfaces. The 
functional relationships between these measures of model complexity and grid 
densities were set heuri stical 1y by varying the grid densities and observing 
the performance of the hidden line routine for several models. For the 
smallest models a grid density of 3 X 3 is used, and for the largest models a 
13 X 13 grid density is used in the JONES-D hidden line routine. For models 
significantly more complex than the largest test models, a higher grid density 
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’ would be required for optimum performance. The lists of vectors and surfaces 

for each of the X-Y cells are then sorted by depth (Z). The end result of the 
spatial sorting are depth sorted cell vector and surface lists. 

♦ 

, 0 EDGE VISIBILITY: A vector (Vj) is sequentially pulled from the global 

vector list. Its location in the spatial map, figure 3 , is determined. The 
list of cell vectors is binary searched to find the depth to limit the search 
for vectors that intersect with Vj. Its intersection with all vectors that 
are ahead of it, are surface edges, and not found to be invisible by a prior 
calculations, are determined. A segmented vector is created from Vj using its 
end points and the calculated points of intersection. Each segment is either 
all visible or invisible. The mid-point of each segment is computed and 
checked against the appropriate cell surface list to ascertain visibility. 

This requires finding the surface cell that contains the vector segment 
mid- point and then performing a binary search to find the depth in the cell 
surface list to limit the search for surfaces that hide the segment midpoint. 
Containment ana depth tests are then performed to ascertain midpoint 
visibility. 

0 RETURN DATA: At the completion of the plot some of the data arrays 

created ■'n JONES-D are passed to the main NPLOT routine to facilitate the 
labeling of the visible grid points and elements. 


Several restrictions on the code, as implemented, should be noted but they in 
general have no impact on the plotting of NASTRAN finite element models. 
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0 A line penetrating a surface results in a visible plot error. This is 
desirable for NASTRAN plotting since this usually indicates a modelling error. 

0 Grid points are required at the points where elements intersect. This 
normally is the case in NASTRAN niodels. 

0 Surface and solid element topologies must be reduced to four node flat 
surfaces. This presents no problems for concnonly used NASTRAN elements. 
Triangles through 20 node brick elements are processed by NPLOT to this format, 


, PERFORMANCE 


The performance of the JONES-D hidden line technique was evaluated several 
different ways: on an absolute basis, in comparison with other plot types, 

and in comparison with other hidden line methods. In order to present an 
accurate picture several things should be noted: 

a. The run tinies (CPU and elapsed) presenteo herein include the time to 
execute the plot function nodule (hidden line, halo line or normal plot), run 
the PLOT 10 module, and paint the picture on a CRT screen. 

b. The CPU times show about a ^5 percent variation due to the work load 


on the VAX. 
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One set of tests consisted of running each of the 14 test models (table 1) to 
generate normal plots, haloed plots and hidden line plots using NPLOT, and 
recording the CPU and elapsed times for each of the plots. To eliminate any 
bias due to direction sensitivity, each run consisted of making three plots in 
orthogonal directions, and averaging the CPU and elapsed time. One measure of 
the performance of JONES-D is presented in figure 4 in the form of a plot of 
the edge vector processing rate (vectors per CPU second) as a function of 
problem size. Each data point represents the average processing rate for one 
of the 14 test models. The overall processing rate appears to be fairly 
linear with an average rate of 103.7 vectors per CPU second. The slight roll 
off at the high end of model size was attributed to the increase page faulting 
generated by these larger models. Recent experience using version 3.0 of the 
VAX operating system and a larger working set size (512 kbytes) increased the 
processing rate for these larger models by about ten percent. 

Presented in table 2 is a comparison of the processing rates for normal (all 
lines visible) plotting, halo plotting, and hidden line plotting. The average 
processing rate for normal plotting was 134.7 vectors per CPU second, for 
hidden line, 103.7, and for haloed line, 55.7. At this point, it should be 
noted that the elenient load section of NPLOT was designed to flag certain 
classes of highly redundant edge vectors that result from solid elements. 

These vectors arc eliminated very early in the JONES-D hidden line routine. 
Thus, these vectors were not counted in calculating the processing speed of 
the hidden line routine. The only two models with solid elements used in the 
effort, MIRROR and hIRRORH, show the effect of this procedure. For example. 
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the CPU tin« to hidden line plot MIRROR was 9.1 seconds. While MIRROR 
nominally consists of 1872 vectors, only 828 vectors need to be fully 
processed by the hidden line routine; this results in a calculated processing 
rate of 91 vectors per CPU second for the hidden line plotting. 

The CPU tine performance of the JONES-D hidden line technique for the 14 test 
models is compared to that for normal and haloed plots in table 3. As was 
expected, halo plotting was the slowest. In most cases, normal plotting was 
slightly faster than hidden line plotting. On an average CPU time basis, 
normal plotting was 16 percent faster than hicden line plotting, and halo 
plotting was 131 percent slower than hidden line plotting. 

To the interactive user, elapsed (wall clock) time can be more significant 
than CPU time. Presented in table 4 is a comparison, based on elapsed time, 
of bidden line plotting, haloed plotting, and normal plotting for the 14 test 
models. In all cases except one, hidden line plotting was the quickest. In 
all cases haloed line plotting was the slowest. On average, the elapsed tinte 
for normal plotting was 15 percent slower than hidden line plotting, and 
haloed line plotting was 198 percent slower than hidden line plotting. 

The JONES-U algorithm was compared to five other hidden line methods. Three 
of the comparisons were based on hands-on use of the alternative algorithms 
and two of the comparisons were based on published performance data. This was 
not meant to be a precise, exhaustive, and detailed comparative study but 
rather to illustrate relatively large differences in performance when dealing 
with the target application, which was NASTRAN finite element model plotting. 



a. The Watkins hidden line/ surface method is used in the MOVIE program. 

A VAX implementation of MOVIE was used for this study. The Watkins method 
does not support line element types so an all surface model, the surface 
elements of the FLFSS model, was used for this comparison. In this form the 
test model contained 857 surfaces. The CPU time for MOVIE was 41.5 seconds 
for generating the hidden line plot; the corresponding time for JOhES-D was 
11.6 CPU seconds. The MOVIE plot contained a few plot errors; however, the 
author's understanding is that the most recent version of MOVIE is fixed in 
this regard but that it runs a little slower. 

b. The SKETCH hidden line routine recently developed by Hedgley 
(reference 12) was obtained and converted to the VAX 11/780 computer. The 
routine as delivered was limited to about 250 polygons, therefore, a 
relatively small model wac selected for testing (the CQUADS of the FSS model, 
183 surfaces). The CPU time for SKETCH was 19.3 seconds for the hidden line 
plot; the corresponding time for JONES-D was 3.5 CPU seconds. This level of 
performance for SKETCH {r^9 polygons per second) seems consistent with the 
data presented in reference 12. In that reference the processing rate for 
SKETCH on a CUC 6500 computer (slightly faster than a VAX 11/780) for finite 
element type models was about 10 polygons per CPU second. Some plot errors 
were observed in the SKETCH generated plot, but this may be due to the fact 
that single precision (32 bit) data formats were used on the VAX; whereas, the 
original implementation used a 60 bit format (CDC 6500). 
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c. The VIEW thermal view factor program developed for NASA/LaRU by 
Profeusor Emery of the University of Washington contains a hidden line 
routine. The NCONES thermal model, figure 5, was used fo>' comparison testing. 
The CPU time for VIEW was 48.8 CPU seconds for the hidden line plot, the 
corresponding time for JONES-D was 8.05 CPU seconds. VIEW was developed on a 
computer (PUP 11/45) with limited memory resources, so by necessity it does a 
lot of FORTRAN I/O which may account for its relatively slow performance. 

VIEW does seem very reliable, and its plots have contained no obvious plot 
errors. 


d. A fast hidcen line method was developed by Bareau (reference 2) for 
the plotting finite element models. We were not able to obtain the code 
for this technique, so no hands-on testing was performed. Based on the data 
presented in the referenced paper it seems quite fast. For example, on a CDC 
6500 computer, (about 1.2 times faster than a VAX 11/780), its performance was 
about 55 edge vectors per CPU second. This levj of performance is 2.5 to 3.0 
times slower than JONES-D. 

e. In a recent paper, (reference 13) Wittram presented a new fast hidden 
line niethod. Again we were not able to obtain the code necessary for hands-on 
testing. For a scene with 1029 faces, the CPU time is given as 20 CPU 
seconds; to this must be added the X-sort time of about 15 CPU seconds for a 
total CPU tine of 35 seconds on a ICL 1906S computer. This corresponds to a 
processing rate of about 30 surfaces per CPU second; the performance of 
JONES-U for a similar size nedel (FLFSS, surface eletrents only) is about 74 
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surfaces per CPU second. Unfortunately, the relative performance of the ICL 
1906S to the VAX 11/780 is somewhat uncertain. Based on Wittram's comment 
that a ICL 1906S is about 10 times slower than an IBM 3081, one could surmise 
that a VAX 11/780 is about 20 percent fastef* than an ICL 1906S. If this is 
true, then Wittram's method is about two times slower than OONES-D. 


COhLL USIOhS 

A high speed hidden line technique has been developed to facilitate the 
plotting of hASTRAN finite element models. Based on testing using many (14) 
different models, the new hidden line algorithm (JONES-D) appeal's to be: 

a. Very Fast: Its speed equals that for normal (all lines visible) 

plotting, and when compared to other existing methods, it appears to be 
substantially faster. 

b. Very Reliable: ho plot errors have been observed using the new 

niethoo to plot MASTRAN models. 

The new 'algorithm nas been made part of the hPLUT NASTRAN plot package and has 
been used at the GSFC by structural analysts for normal production tasks. It 
is interesting to note that analysts when given a choice, seldom use normal, 
all lines visible, plotting. 
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MODEL 

GRID 

POINTS 

LINE 

ELEMENTS 

SURFACE 

ELEMENTS 

SOLID 

ELEMENTS 

JPL 

55 

9 

67 

0 

FSS 

167 

352 

279 

0 

CDM 

210 

180 

191 

0 

MAYPOLE 

231 

542 

96 

0 

OSSFF 

315 

462 

130 

0 

CUBE8D 

408 

382 

342 

0 

BBXRT 

535 

0 

576 

0 

FLFSS 

604 

1235 

857 

0 

MiMS 

678 

678 

519 

0 

DHRS 

703 

717 

779 

0 

MIRROR 

540 

0 

180 

720 

FEM 

1268 

1146 

1012 

0 

LSDBUCK 

1502 

1683 

1382 

0 

MIRRORH 

1183 

0 

0 

864 


NASTRAN Model Data 


Table 1 


















SPEED VECTOR/CPU SEC. 


MODEL 

VECTORS 

NORMAL 

HALO 

HIDDEN 

JPL 

126 

131.3 

68.1 

97.6 

FSS 

402 

128.8 

58.6 

106.1 

CDH 

472 

135.2 

81.4 

131.5 

MAYPOLE 

518 

136.3 

15.4 

78.4 

OSSFF 

525 

127.4 

65.1 

147.5 

CUBE8D 

894 

134.4 

59.8 

124.7 

BBXRT 

1084 

139.0 

60.5 

97.7 

FLFSS 

1375 

133.5 

55.1 

98.4 

MJ>1S 

1430 

137.5 

62.6 

100.3 

DHRS 

1561 

136.8 

53.9 

102 

MIRROR 

1872/828* 

133.1 

31.7 

91.0 

FEM 

2413 

136.6 

46.1 

88.1 

LSDbUCK 

3056 

141.1 

69.5 

82.7 

mirrokh 

3198/1152* 

135.2 

52.1 

106.5 


* Nuniier of vector for hidden line processing 


Vector Processing Speed 


Table 1 


hODEL 

VECTORS 

NORMAL 

HALO 

HIDDEN 

JPL 

126 

.96 

1.85 

1.29 

FSS 

402 

3.12 

6.86 

3.79 

CUH 

472 

3.49 

5.80 

3.59 

MAYPOLE 

518 

3.80 

33.59 

6.61 

OSSFF 

525 

4.12 

8.07 

3.56 

C0BE8D 

894 

6.65 

14.96 

7.17 

BBXRT 

1084 

7.80 

17.91 

11.10 

FLFSS 

1375 

10.30 

24.93 

13.97 

MMS 

1430 

10.40 

22.86 

14.25 

DHRS 

1561 

11.41 

28.96 

15.31 

MIRROR 

1872/828* 

14.06 

58.96 

9.10 

FEM 

2413 

17.66 

52.33 

27.39 

LSDBUCK 

3056 

21.66 

43.99 

36.96 

MIRRORH 

3198/1152* 

23.65 

61.35 

10.82 


* Number of vectors for hidden line processing 


CPU Time Comparison 
Hidden vs. Halo vs. Normal 
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hODEL 

VECTORS 

t 

OPL 

126 

FSS 

1 

402 

CDH 

472 

MAYPOLE 

518 

OSSFr 

525 

C0BE8D 

894 

1 BBXRT 

1084 

! FLFSS 

1375 

j| MFiS 

1430 

DHRS 

j 

1561 

1 MIRROR 

1872/828* 

j FEM 

1 

2413 

1 

1 LSDBUCK 

3056 

i MIRRORH 

3198/1152* 


ELAPSED TIKE (SEC.) 


NORMAL 

HALO 

HIDDEN 

3.6 

5.3 

3.2 

11.2 

23.9 

8.0 

12.5 

19.7 

7.8 

13.3 

84.6 

17.4 

13.9 

22.6 

8.9 

25.2 

37.5 

14.8 

27.1 

52.1 

23.5 

35.8 

65.4 

25.4 

34.2 

59.9 

27.5 

40.0 

70.2 

26.9 

47.9 

214.8 

19.0 

61.0 

124.9 

60.0 

74.2 

86.6 

68.7 

85.5 

131.6 

23.6 


i 

\ 
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* Number of vectors for hidden line processing 


Elapsed Time Comparison 
Hidden vs. Halo vs. Normal 



Table 4 
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FIGURE 2 
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Presented in this appendix are the hidden line plots for the fourteen test 
models (figures lA through 14A) generated by use of the JONES-D hidden line 
routine and the NPLOT NASTRAN plot package. Noted on each plot is the CPU 
time required to generate it. Illustrated in figure 15A is the effect of 
element labeling and grid labeling on a hidden line plot. On complex plots, 
the zoom feature is essential to declutter the displayed image;- presented in 
figure 16A is such a zoom view of the image shown in figure 15A. 
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SOURCE LISTING 
dONES-D HIDDEN LINE ROUTINE 
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Presented in this appendix are the source listings of: 

a. The JONES-0 hidden line routine, pages B2-B11. 

b. The code used to generate the SORTP array in the main NPLOT routine, 
page B12. 

c. The SORTLEN subroutine called by the hidden line routine, page B12. 

d. The VEC subroutine called by the hidden line routine to plot absolute 
vectors, pages B13-B14. 
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HIDDEN LINE MODULE: HIDDEN 


This module uses the JQNES-D hidden line technique 
developed by Gary Jones# NASA/GSFC# Mail Code 731# 
Greenbelt, MD# 20771 <Tel: 301 344-7166) 

This module is designed to function with the NPLOT 
Nastran plot package. 


VERSION DATE: April 20# 1902 


PARAMETER DESCRIPTION: 

PGRID — Coordinate locations of grid points afte^ rotations 
and perspective transformation. 

VECLIS - Global vector list# pointers to PGRID. 


SURLIS - Global surface list# pointers^ to PGRID. 

NVEC - Number of vectors in the global vector list# VECLIS. 
NSUR - Number of surfaces in the global surface list. SURLIS 


X1#Y1#Z1 “ View parameters: 2# 1#3 — XY view 

3. 2# 1 = Y2 view 
1# 3# 2 «= XZ view 


VISLIB - Vis-ibility 
visibility; 
subsequent 


table for grid points. A value of 2 indicates 
passed to the main NPLOT routine to faciltate 
labelling of visible grid points and elements. 


XMAX# XMIN# YMAX# YMIN - Horizontal and vertical max and min values 


SORTP — Shell sort parameters# set up in NPLOT main. 

SOLIDVV - Counter used to reject highly redundant solid 

edges. Set up is accomplished at time of element load 

in NPLOT. 


NSVEC - Number of vectors in the global vector list that are 
edges of surface or solid elements. 


XYS# NUMS, SURF# IBUCK5# XDELT5. YDELTS - Surfac 
routine; passed to the main NPLOT 
subsequent labelling of visible e 


e data created in this 
routine to facilitate 
Icmen ts. 


HP - Flag to indicate if the plot file option is selected. 

a file that captures the on screen plot data is Senerated. 

The HP4 software uses this file to create plots on an HP7580A 
pen plotter. 


VECTYPE - Counter used to indicate that vectors result solely 

from line element types. If VECTYPE< I ) *=0 the edge vector 
is not a surface or solid element edge. 
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SUBROUTINE H1DDEN(PGRID. VECLIS. SURLIS. NVEC. NSUR. XI. Yl. Zl. VISLIS 
&. YMAX. YMIN. XMAX. XHIN. 50RTP. SOLIDW. NSVEC 
&. XYS. NUMS. SURF. 1 DUCKS. XDELTS. YDELTS. HP. VECTYPE) 

INTEGER»2 VECL1S<9000. 2) . SURLIS (9000. 4 ) , VISLIS<9000) 

&. XYV(600. 13. 13). NUt1V< 13. 13). VP<9000. 4). XYS(600. 13. 13). NUHS(13. 13) 
&. JUMP. UMAX. J3. J4. SORTP( 8). VECTYPE (9000) 

&. XI. Yl. Zl. SOLIDW (9000). VECVP( 9000) 

REAL*4 PGRID(9000. 3) . VECTOR ( 1 5. 9000) . SURF (26. 9000) 

&. XMAX. XMIN. YMAX. YMIN. LINE( 1000. 4) 

LIMNVEC-600 

LIMNSUR=600 

C * 

919 FORMAT OX. 2E1 3. 6) 

c ♦ 

C ♦ PREP WORK FOR HIDDEN LINE ♦^M^*♦♦*♦♦*♦♦♦♦♦«*#**#«♦^M^^^***♦***♦*♦**♦ 
C * 

C « CALCULATE GRID DENSITIES FOR X-Y BUCKET SORTS 
IF(NSVEC. LT. 200) THEN 
IBUCKV=3 
GOTO 3 

ELSE IF (NSVEC. LT. 600) THEN 

IBUCKV«6 

GOTO 3 

ELSE IF (NSVEC. LT. 1500) THEN 

IBUCKV=»8 

GOTO 3 

ELSE IF (NSVEC. LT. 2400) THEN 

IBUCKV*10 

GOTO 3 

ELSE 

IDUCKV«13 

ENDIF 

3 XYVFAC = IDUCKV-, 01 

IF (NSUR. LT. 100) THEN 

IBUCKS=3 

GOTO 4 

ELSE IF (NSUR. LT. 600) THEN 

IBUCKS=G 

GOTO 4 

ELSE IF (NSUR. LT. 1200) THEN 

IBUCKS^IO 

GOTO 4 

ELSE 

IBUCKS«13 

ENDIF 

4 XYSFAC«IDUCKS-. 01 
XDIF~(XMAX-XMIN> 

IF(XDIF. EQ. 0. 0) THEN 
XDIF = . 01 

ENDIF 

YDIF = (YMAX’^YMIN) 

IF(YDIF. EQ. 0. 0) THEN 

YDIF» . 01 

ENDIF 

XDELTV=XYVFAC/XDIF 

XDELTS=XYSFAC/XDIF 

YDELTV=XYVFAC/YDIF 

YDELTS=XYSFAC/YDIF 
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DO 5 I»l, lOUCKV 
DO 5 J=l, IBUCKV 

5 NUMV(J. I)-0 
DO 6 I»l. IBUCKS 
DO 6 J=l. IBUCKS 

6 NUMS(J. I)-0 

C * 

C * BEGIN VECTORILINE) PREP ••»«*«*««««•«•* 

C * * 

DO 140 NPT=1,NVEC 
IF(SOLIDW<NPT). OT. 23) GOTO 140 
VECVPINPT)*! 

Nl»VECLIS(NPT, 1) 

N2«»VECLIS<NPT, 2) 

VECT0R(1,NPT)=PGRID<N1. XI) 

VECT0R(2. NPT)«=PGRID(NI. Y1 > 

VECTORO. NPT)«=PGRID(N1. Z1 ) 

VECTaR(4, NPT)=PGR 1 D(N 2 . XU 
VECT0R(5, NPT)=PGRID(N2. Y1 ) 

VECT0R(6, NPT)“PGR1D(N2. Z1 ) 

VECTOR ( 7. NPT) “VECTOR (5. NPT) -VECTOR (2. NPT) 

VECTORIQ, NPT)=VECTOR( 1, NPT)-VECT0R(4. NPT) 

VECTOR (9. NPT ) “-VECTOR (7t NPT ) *VECTaR < 4. NPT) -VECTOR < 0 . NPT) 
8«*VECT0R(5. NPT) 

VECTOR( 10. NPT) “MAX I VECTOR! 1. NPT), VECTOR (4. NPT) ) 

VECTOR (11, NPT) “MAX (VECTOR (2, NPT) . VECTOR (5. NPT) ) 

VECTOR <12. NPT) “MAX (VECTOR (3. NPT ), VECTOR (6. NPT) ) 

VECTOR! 13, NPT)=MIN(VECTOR( 1. NPT), VECT0R(4, NPT) ) 

VECTOR! 14. NPT)-MIN!VECT0R(2. NPT ). VECTOR ! 5. NPT) ) 

VECTOR (15. NPT) “MIN! VECTOR (3, NPT ), VECTOR ( 6, NPT) ) 

C * 

C * BUCKET SORT !X.Y) ON VECTORS TO PRODUCE VECTOR MAP !XYV) 
ISTX=1 + (VECT0R(I3. NPT)-XMIN)*XDELTV 
ISPX“lf!VECTOR! 10. NPT)-XM1N)*XDELTV 
ISTY“l-t-( VECTOR (14, NPT ) -YMI N > *YDELTV 
ISPY=1 ♦! VECTOR! 11, NPT) -YMIN)*YDELTV 

C * 

IF!VECTYPE!NPT). CE. 1 ) THEN 
DO 137 J-13TX. ISPX 
DO 137 1“ISTY. ISPY 
NUMV < J, 1 ) “NUMV < U, I ) ••■ 1 

137 XYV(NUMV(U, I ). J. 1 )=NPT 
ENDIF 

C ♦ 

138 CONTINUE 
VP(NPT. 1 )=ISTX 
VP!NPT, 2>“ISPX 
VP(NPT. 3)“ISTY 
VP(NPT. 4)=ISPY 

140 CONTINUE 

C • 

C * SHELL SORT BY DEPTH (2) OF VECTOR MAP (XYV) 

DO 145 IX“1, IBUCKV 
DO 145 IY-1. IBUCKV 
JLEN=NUMV( IX, lY) 

IF<JLEN LE 1) GOTO 145 
CALL SORTLENIULEN, UCT) 

DO 144 JINDEX-JCT, 1 , -1 
JUMP“SORTP (UINDEX ) 

UMAX “JLEN- JUMP 

146 IFLIP=-0 
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DO M3 M-l.JMAX 
N-M+JUMP 
J3«»XYV(li. IX, IY> 
j 4«XYV(N. IX. lY) 

IF<VECTDR( 15. J3). CT. VECTOR! 15. J4 ) ) THEN 

XYV(M. IX. IY)»J4 

XYV(N. IX. IY)»J3 

IFLIP-1 

END IF 

143 CONTINUE 

IF<IFLIP. EQ. 1) GOTO 146 

144 CONTINUE 

145 CONTINUE 

C • 

C ♦ CHECK FOR OVERFLOW IN VECTOR MAP, XYV 
DO 149 I-l, IBUCKV 
DO 149 U“1,IBUCKV 
IF(NUMV< J, I ). GT. LIMNVEC) THEN 

TYPE 147. LIMNVEC 

147 FORMAT!/. IX. 'OVERFLOW IN VECTOR MAP XYV'.9X. 'LIMIT 

GOTO 800 
END IF 

149 CONTINUE 

C * 

C * BEGIN SURFACE PREP ************** 

C * 

DO 150 NPT=1.NSUR 
N1»SURLIS!NPT. 1) 

N2“SURLIS(NPT. 2) 

N3“SURLIS!NPT. 3) 

N4=SURLIS!NPT. 4) 

XP1«PCRID<N1. XI ) 

XP2-PGRID!N2. XI » 

XP3*PCRID!N3. XI ) 

XP4-P0RID(N4. XI ) 

YPl-PCRIDCNl. Y1 ) 

YP2=PGR1D!N2. Y1 ) 

YP3-PGRID(N3. Y1 ) 

YP4-PGR1D! N4, Y1 ) 

ZP1«PCRID(NI. 21 ) 

2P2-P0R ID!N2. 21 > 

ZP3*PGRID!N3, 21 > 

ZP4«PCRID(N4. 21 > 

SURF!1.NPT)“XP2-XP1 
0URF!2. NPT)-XP3-XP2 
SURF!3. NPT>«XP4-XP3 
SURF! 4, NPT)“XP4-XP1 
SURF! 5, NPT)»XP1-XP3 
SURF!6. NPT)»YP1-YP2 
SURF ! 7, NPT ) -VP2-YP3 
SURF!8. NPT)=YP3-YP4 
SURF! 9, NPT)=YP1-YP4 
SURF < 1 0, NPT ) =YP3-YP 1 
SURF! 1 1. NrT)=XP2*YP3-XP3*YP2 
SURF! 12, NPT)=XP3*YP1-XP1*YP3 
SURF! 13. NPT)-XPl*YP2-XP2*YPl 
SURF! 14, NPT)=XP3»YP4-XP4*YP3 
SURF! 15. NPT)=XP1*YP4-XP4*YP1 
SURF! 16. NPT)“MAX!XP1, XP2. XP3. XP4) 
me / t -7 fcjOT ^wMTMfXPl. )CP2» XP3» XP4) 


14) 
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B6 SURF( 19, NPT)«MIN<YP1, YP2, YP3, YP4) 

BURF(20, NPT)«MAX(2PJ, ZP2, ZP3. 2P4) 

SURF(21#NPT)«MIN(ZP1, ZP2, ZP3. ZP4) 

C * 

C ♦ BUCKET SORT <X,Y) ON SURFACES TO PRODUCE SURFACE MAP (XYS) 
ISTX»l + (SURF<17. NPT)-XMIN)*XDELTS 
ISPX=1-KSURF( 16, NPT)-XMIN)*XDELTS 
ISTYo 1 ( SURF ( 1 9, NPT ) - YMIN ) * YDELTS 
ISPY*=1 + (SURF<18. NPT)-YMIN)*YDELTS 
335 DO 337 J=ISTX. ISPX 
DO 337 I«ISTY, ISPY 
NUMS(Ji I)«NUMS(J* 

337 XYS(NUMS(J, I), J, 1)«NPT 

C * 

C SURFACE AREA CALCULATION 

AREA=ADS(SURF< 11. NPT)-^SURF( 12» NPT ) ♦SURF ( 13t NPT) ) 
«<^ADS(SURF( 12, NPT) + XPl-»YP4-t-XP4*SURF( 10, NPT )-XP3# YP4 ) 

IF(AREA. EQ. 0. 0) THEN 
SURF(22, NPT)»-1. 0 
ELSE 

SURF ( 22, NPT ) -■ 1 /AREA 
ENDIF 

C * 

C ♦ SET SURFACE ELEMENT TO ZERO 
SURF(25. NPT)*=0. 0 

150 CONTINUE 

C ♦ 

C * SHELL SORT BY DEPTH <Z) OF SURFACE MAP (XYS) 

DO 450 IX«=1, IBUCKS 
DO 450 IY=1. IBUCKS 
JLEN=NUMS( IX, lY) 

IF(JLEN. LE. n GOTO 450 
CALL SORTLEN( JLEN, JCT) 

DO 435 JINDEX-JCT, 1,-1 
UUMP«SORTP ( JINDEX ) 

JMAX«JLEN-JUMP 
446 IFLIP-0 

DO 400 M«1,JMAX 
N»M^JUMP 
J3«XYS(M, IX, lY) 

J4«XYG(N, IX, lY) 

IF(5URF(21, J3). GT. SURF(21, J4) ) THEN 

XYS(M, IX, IY)«J4 

XYS(N, IX, IY)*=J3 

IFLIP«1 

ENDIF 

400 CONTINUE 

IF( IFLIP. EQ, 1 ) GOTO 446 
435 CONTINUE 
450 CONTINUE 

C * 

C * CHECK FOR OVERFLOW IN SURFACE MAP. XYS 
DO 156 1 = 1, IBUCKS 
DO 156 J=l, IBUCKS 
IF(NUM5( J, I ). GT. LIMNSUR) THEN 
TYPE 157,LIMNSUR 

157 FORMAT(/, IX, 'OVERFLOW IN SURFACE MAP XYS ' , 9X , 'L I M I T =',I4) 
GO TO 800 
ENDIF 

156 CONTINUE 
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607 

C * 

C ♦ 

615 

608 

C * 
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DO 800 IbI.NVEC 
1F(S0LIDW<I). CT. 23) GOTO 800 
FLIP*0 
NSEG-1 

ICPA»VECLIS(I. 1) 

IGPB-VECLIS(I,2) 

ICHK=1 

IF<ABS(VECT0R<7, I) ). GT. ABS(VECT0R(8. I > > ) ICHK»2 
IF<VECT0R(7, I). EQ. 0. 0. AND. VECT0R(8, I ). EQ. O. 0) GOTO 800 

THE ARRAY LINE(I.l-3) HOLDS THE VECTOR END POINTS 
AND POINTS OF INTERSECTION. LINE(I>4) INDICATES 
SEGMENT VISIBILITY! O » HIDDEN. J - VISIBLE 
LINEd. 1)«VECT0R(1, I) 

LINE(1,2)=VECT0R(2. I) 

LINEd, 3)-VECT0R(3. I) 

LINEd, 4) = 1. 

XMINI=VECT0R(13. I) 

YMINI-VECTQRd4. I) 

XMAXI-VECTORdO, I) 

YMAXI=VECTORdl, I) 

ZMAXI»VECT0Rd2, I) 


87 


LINE INTERSECTION CALCULATION ••»•««»«*«•«»«»«• 

DO 610 IX=VP( I, 1 ), VPd, 2) 

DO 610 IY-=VP<I,3), VPd,4) 

BINARY SEARCH TO OBTAIN SEARCH DEPTH IN VECTOR MAP 
LEN-NUMV(TX. lY) 

IF(LEN. LT. 3) GOTO 615 

MAXX-LEN-l 

LOW-0 

LAST=0 

JPC-LEN/2 

DO 607 WHILE (LAST. NE. JPC ) 

LAST-JPC 

J3»XYV( JPC, IX. lY) 

IFIVECTORdS, J3). LT. ZMAXI ) THEN 
J3=XYV(JPC+1. IX. lY) 

IF(VECT0Rd5, J3). GE. ZMAXI) GOTO 608 

LOW-JPC 

ELSE 

MAXX-JPC 

ENDIF 

JPC=(LQW+MAXX+1 )/2 
CONTINUE 


DETERMIN'l LINE INTERSECTIONS 
1F(JPC. EQ. 1) GOTO 610 
JPC»LEN 

DO 609 JP=1,JPC 
J=>XYV( JP, IX. IY> 


VECVP<J)=0 IF VECTOR J HAS BEEN FOUND TO BE INVISIBLE 
IF(VECVP(J). EQ. 0) GOTO 609 

IF ( VECTOR { 13, J>. LE. XMAXI. AND. VECTOR ( 1 0. J > . CE. XMINl . AND. 

fcVECTOR { 14, J). LE. YMAX I. AND. VECTOR ( 1 1 , J) . GE. YMINI ) THEN 
IF( VECLIS( J, 1 > . EQ. VECLIS< 1 , 1 ) . OR. VECL I S ( J, 1 ) EQ VECLIS<I,2) 
^^°^^'^^CLIS(J. 2). EQ. VECLIS( I. 2). OR. VECLIS( J. 2). EQ. VECL1S( I. 1 ) ) 

DEN0M=VECT0R(7, I )*VECTOR(0. J) -VECTOR ( 7, J ) *VECTOR <0, I » 



OniGlNAl PAGH 

B8 OF POOR quality 

1F(DEN0M. NE. 0. 0) THEN 
IFdCHK. EQ. n THEN 

XI NT- (VECTOR <8. I >*VECT0R<9. J)-VECT0R<8. J)*VECT0R<9, I) >/DENOH 
IFCXINT. OT. XMAXI. OR. XINT. LT. XMINl > GOTO 609 
VX1-<VECT0R<1. I)-XINT)/VECT0R<8. I> 

YINT-VECT0R<2. I >-«-VXl*VECT0R<7, I) 

ZINT-VECT0R(3. !)♦( VECTOR (6. I) -VECTOR (3, I)>*VXI 
ELSE 

Y I NT- ( VECTOR (9.1) ^VECTOR ( 7. J ) - VEC TOR ( 7, I ) * VECTOR ( 9. J ) ) /OENOM 
IF(YINT. OT. YMAXI. OR. YINT. LT. YMINI ) GOTO 609 
VYI»( VECTOR (2. I >-YINT> /VECTOR ( 7. I ) 

XINT»VECT0R(1. I >^VYI*VECT0R(8. I ) 

ZINT-VECT0R(3. I)-KVECT0R(3. I)-VECT0R(6. I))*VYI 
END IF 

IF(ABS(VECT0R(8, J) ). OT. ABS(VECTOR ( 7. J) ) ) THEN 
IF(XINT. OT. VECTOR( 10, J). OR. XINT. LT. VECTOR < 13. J) ) THEN 
XITF-. 001*(ABS(XINT) ) 

XIT-XINT-XITF 

IF( X IT. OT. VECTOR (10. J) ) GOTO 609 
XIT-XINT+XITF 

IF(XIT. LT. VECTOR! 13. J) > GOTO 609 
ENDIF 

ZJNT-VECT0R(3, J)-*-(VECT0R(3. J)-VECT0R(6. J) )* 

«c(XlNT-VECTOR(l. J) )/VECT0R(8, J) 

ELSE IF (VECT0R(7. J). NE. O. O) THEN 

IF(YINT. CT. VECTOR! 11. J). OR YINT. LT. VECTOR! 14. J) ) THEN 
YITF-. 001*(ABS(YINT) ) 

YIT-YINT-YITF 

IF(Y1T. GT. VECTOR! 11. J)) GOTO 609 
YIT-YINT+YITF 

IF!YIT. LT. VECTOR! 14. J) ) GOTO 609 
ENDIF 

ZJNT=VECT0R!3, J) + !VECT0R(3, J)-VECT0R!6, J) »♦ 

«<!YINT-VECT0R!2, J) )/( -VECTOR ( 7, J) ) 

ELSE 

GOTO 609 
ENDIF 

C * 

C * DEPTH AT INTERSECTION TEST 
IF!ZVNT LT. ZINT) THEN 
NSEG-NSEGfl 
LINE!NS£G. 1 )=X1NT 
LINE!NSEC. 2>“YINT 
LINE(NSEG. 3)=ZINT 
LINE!NSEG. 4)=1. 

ENDIF 

C ♦ 

ENDIF 

ENDIF 

609 CONTINUE 

610 CONTINUE 

C * 

C ♦ INSERT END POINT 
N5EC»N5EG^1 

LI NEC NSEC, 1 )=VECT0R(4. I ) 

LINECNSEC, 2)»VECT0R<5, I> 

LINECNSEG, 3)*VECT0R(6, I) 

LINECNSEC. 4)=1. 

C ♦ 

C » SORT INTERSECTION LIST FROM MIN TO MAX 
IFCNSEG LE. 3) GOTO 600 
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1F<LINE( X. ICHK). CT. LINE(NSEC. ICHK) > FLIP-X 
IFCNSEO. LE. 5) THEN 
JCT-X 
GOTO 62 X 

ELSE IF (NSEC. LE. X3» THEN 

JCT-2 

GOTO 62X 

EL^^IF (NSEC. LE. 29) THEN 

JCT-3 

GOTO 62X 

ELSE IF (NSEC. LE. 61) THEN 

JCT-4 

GOTO 62 X 

ELSE IF(NSEG.LE. X25) THEN 

JCT-5 

GOTO 62 X 

ELSE 

JCT-6 

ENDIF 

62X CONTINUE 

DO 620 UINDEX-JCT, X. -1 
UUMP-SORTP( JINDEX) 

JMAX-NSEC-UUMP 
626 IFLIP-0 

DO 623 MaX.JMAX 
N=MfJUMP 

IF(LINE(M, ICHK). OT. LINE(N. ICHK) ) THEN 
DO 62G J3-X. 4 
TEMP»LINE(M. J3) 

LINE(M. U3)«LINE(N. J3) 

LINE(N. J3)-TEMP 
628 CONTINUE 
IFLIP-X 
ENDIF 

623 CONTINUE 

IF( IFLIP. EQ X ) GOTO 626 
620 CONTINUE 
600 CONTINUE 

c * 

C * COMPUTE VISIBILITY OF LINE SEGMENTS 
DO 790 J»2, NSEG 
XMID-<LINE(U-X, X INLINE (U. 1 ) > /2 
YMID-(LINE(U-X. 2)>LINE(U. 2) )/2 
2MID-(LINE(J-1.3)-»L1NE(U. 3) )/2 
I X = l-*- ( XMID-XMIN)«XDELTS 
IY“t + < YMID-YMIN)*YDELTS 
c • 

C • DINARY SEARCH TO OBTAIN SEARCH DEPTH IN SURFACE MAP 
LEN*=NUMS( IX. :Y) 

IF(LEN LT. 3) GOTO 715 

MAXX«LEN-1 

L0W«0 

LAST-0 

JPC-LEN/2 

DO 751 WHILE <LAST. NE. JPC) 

LAST-JPC 

J3-XYSCJPC. IX. lY) 

IF(SURF(2l. J3). LT. 2MID) THEN 
J3 = XYSl JPC'^l. IX. lY) 

IF(SURF(21. U3). GE, 2MID) GOTO 753 
LOW-JPC 
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ELSE 

MAXX-JPC 

ENDIF 

JPC“<L0W+MAXX+n/2 
7S1 CONTINUE 

IF( JPC. EQ. 1 ) GOTO 790 
715 UPC»LEN 

C * 

C ♦ LINE SEGMENT MID-POINT VISIBILITY TESTS 

753 DO 709 JP=1.JPC 
U3=XYS( JP. IX. IY> 

IF(SURF( 16. J3). GE. XMID. AND. SURF( 17, J3). LE. XMID. AND. 

&SURF< 18, J3). GE. YMID. AND. SURF( 19, J3). LE. YMID) THEN 

C * 

C * CONTAINMENT TEST 

IF(SURF<22. J3). GT. O. O) THEN 

A1=ABS(XMID*SURF(7. J3)+YMID*SURF(2. J3)+SURF< 11. J3) ) 
A2=ABS(XMID*SURF<a, J3)+YMID*SURF(3, J3)+SURF(14, J3) » 
A3=ADS(XMID*SURF<6. J3>+YMID*SURF( 1. U3)+SURF( 13, J3> ) 
A4=ADS(XMID*SURF<9. J3)+YMID*SURF<4. J3)+SURF(15. J3) ) 
ARATI0=*(Al-*-A2+A3-*-A4)*SURF(22. J3> 

IFIARATIO. LT. (1. 001 ) ) THEN 

C * 

C * INITIAL DEPTH TEST 

IFIZMID. GT. SURF(20. J3) > GOTO 780 

C * 

C * IS MID-POINT ON THE SURFACE BOUNDARY? 

IF( (A1*SURF(22. J3) ). LT. (. 0005) ) GOTO 789 
IF( (A2»SURF(22. J3) ). LT. (. 0005) ) GOTO 789 
IF( <A3*SURF(22. J3) ). LT. (. 0005) ) GOTO 789 
IF(SURLIS(J3. 1). EQ. SURLIS(J3, 4) ) GOTO 754 
IF< (A4*SURF<22. J3) ). LT. (. 0005) ) GOTO 789 

754 CONTINUE 

IF( IGPA. EQ. SURLIS(U3. 1 ). OR. IGPA. EQ. SURLISI J3. 2). OR. 

«<IGPA. EQ SURLIS(J3. 3). OR. IGPA. EQ. SURLIS< J3, 4) ) THEN 
IF( IGPB EQ. SURLIS( J3. 1 ). OR. IGPB. EQ. SURLIS< J3. 2). OR. 

MGPB. EQ. SURLISI J3. 3). OR. IGPB. EQ. SURLIS(U3. 4) ) GOTO 709 
ENDIF 

C ♦ 

C * IF SURFACE DATA NOT ALREADY COMPUTED THEN CALCULATE IT 
IF<SURF(25. U3). EQ. O. O) THEM 

SURF<25. J3)»SURF( 1 1. J3)+SURF< 12. J3 ) +5URF < 13. U3) 

IF(SURF(25. J3). EQ. O. O) GO TO 709 
ZP1=PGRID(SURLIS( J3, 1 ), Z1 ) 

2P2=PGRID(SURLIS(J3, 2), Z1 > 

ZP3=PGRID(SURLIS( J3. 3). ZD 

SURF(23. U3)=ZP1*SURF(7. 03 ) +ZP2*SURF ( 10. J3) + ZP3*SURF(6. 03) 
SURF(24. 03)=ZP1«SURF(2. 03 ) +ZP2*SURF < 5, 03 > +2P3*SURF< 1 , 03) 

SURF(26, 03)»ZP1*SURF< 11. 03 ) + ZP2*SURF < 12, 03 ) ■«-ZP3*SURF ( 13, 03) 
ENDIF 

C * 

C ♦ DEPTH TEST: IF MID-POINT IS BEHIND SURFACE. THEN VECTOR 

C * SEGMENT IS NOT VISIBLE 

ZSUR=(SURF(26. 03) + XMID*SURF<23. 03 )+YMID*SURF<24, 03) ) /SURF(25. 03) 
IFCZMID GT. ZSUR) THEN 
780 LINE(0. 4)=0. 

GO TO 790 

C * 

ENDIF 

ENDIF 
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ENDIF 

789 CONTINUE 

790 CONTINUE 

C * 

C ♦ 

C * DRAW VISIBLE LINE SEGMENTS 

C * 

VECVP(I)«0 
DO 799 J-2, NSEC 
IF(LINE( Ji 4). NE. 0. ) THEN 

C ♦ VEC - A ROUTINE TO DRAW AN ABSOLUTE VECTOR 

CALL VEC(LINE(J-1, l),LINE(J-l. 2)»LINE(J. 1)#LINE(J. 2) ) 
VECVPCI)»1 

C ♦ IF HP=1 WRITE A PLOT FILE FOR HP75B0A PEN PLOTTER 
IF(HP. EQ. 1> THEN 

WRITEC19.919) LINE( J-1 . 1) . LINE( J-1 . 2) 

WRITE! 19, 919) LINE(J# 1),LINE(J. 2) 

ENDIF 

ENDIF 

799 CONTINUE 
C ♦ 

c ♦ 

c ♦ CHECK AND RECORD VISIBILITY OF END POINTS 

C ♦ 

IF<NSEC. LT. 3. OR. FLIP. EQ. O) THEN 
IF(LINE<2, 4). EQ. 1. ) THEN 
VISLIS(VECLIS(I, 1))=2 
ENDIF 

IF<LINE<NSEG, 4). EQ, 1. ) THEN 
VISLIS(VECLIS(I, 2) )»2 
ENDIF 
ELSE 

IF(LINE(2, 4). EQ. 1. ) THEN 
VISLIS(VECLIS( J, 2) )-2 
ENDIF 

IF(LINE(NSEG, 4). EQ. 1. ) THEN 
VISLIS(VECLIS<I, 1) )*2 
ENDIF 
ENDIF 

C ♦ 

800 CONTINUE 


RETURN 

END 
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CALCULATION OF SORTP . 


C * 

C SORT CRSORT - SHELL METHOD 

C * 

SORTP(l)«l 
DO 431 I « 2, a 
SORTP( I )k1+SORTP( I-X )*2 
451 CONTINUE 


SOR7LEN SUBROUTINE: 


C ♦ 

C * 

C ♦ CALCULATES SHELL SORT LENGTH PARAMETER 

SUBROUTINE SDRTLEN(NLEN, JCT) 

IF(NLEN. LE. 5) THEN 

JCT-1 

GOTO 100 

END IF 

IF(NLEN, LE. 13) THEN 

JCT>t2 

CCTO 100 

END IF 

IF (NLEN. LE. 29) THEN 

JCT«3 

GOTO 100 

ENDIF 

IFINLEN. LE. 61 ) THEN 

JCT«4 

GOTO 100 

ENDIF 

IF (NLEN. LE. 125) THEN 

JCT=5 

GOTO 100 

ENDIF 

IF (NLEN. LE. 253) THEN 

UCT=6 

GOTO 100 

ENDIF 

IF (NLEN. LE. 510) THEN 

JCT=7 

GOTO 100 

ENDIF 

JCT=8 

100 CONTINUE 
RETURN 
END 


c « 
c * 
c * 
c * 
c • 
c * 
c * 
c * 
c ♦ 
c * 

SUBROUTINE VEC(CX1/ CYl/ CX2» 0Y2) 

DIMENSION BUFIN<4), BFOUT(4) 

COMMON /TKTRNX/ TMINUX, TMINVY. TMAXVX, TMAXVY, TREALX, TREALY. 

+ TIMACX. TIMACY, TRCOSF, TRSINF, TRSCAL. TRFACX, TRFACY. 

+ TRPARi, TRPAR2, TRPAR3. TRPAR4, TRPAR5, TRPAR6, KMOFLC<2>, 

+ KGNMOD. KPAOU, KACHAR. KOBLEN. KTRAIL. KLEVEL. KPAD2. 

+ KBAUDR, KCNFLC, KGRAFL. KHOMEY, KKMODE. KHORSZ/ KVERSZ. KTBLSZ. 

+ KSIZEF. KLMRGN, KRMHCN. KFACTR/ KTERM, KLINE, KZAXIS, KBEAMX, KBEAMY. 
+ KMOVEF. KPCHAR(5).KDASHT, KMINSX.KMINSY, KMAXSX. KMAXSY. KEYCON, 

* KINLFT, KOTLFT.KUNIT 

EQUIVALENCE (BFC5UT(1 ), CXS). (BF0UT(2), CYS), (PFOUTI.l), CXE). 
+<BFOUT(4>.CYE) 

C * 

C * THIS SECTION REQUIRED FOR GENISCO G-IOOO TERMINAL 

C * 

CALL BUFFPK<0. ITEMP) 

C * 

C * SET UP OPERATIONS 
C ♦ 

IF . IRAFL. EQ. O) THEN 

TREA. -TMINVX-^(FL0AT(KBEAMX-KMINSX)/TRFACX) 

TREAl . ' •• MINVY+ (FLOAT (KDEAMY-KMINSY > /TRFACY ) 

TIMAi'/ . EALX 
TIMAwV^TREALY 
KGRAFL-1 
ENDIF 

CAU- DUFFPK(1.29» 

KMOVEF-1 

C ♦ 

C • PROCESS START POINT (GX1,CY1) 

C * 

IF(OXl. QT. TMAXVX. OR. GXl. LT. TMINVX. OR. 

■••OYl. OT. TMAXVY. OR. GYl. LT. TMIfA/Y) THEN 
CALL XYCNVT( KBEAMX, KBEAMY) 

ELSE 

DX=GX1-TMINVX 

DY-OYl-TMINVY 

IX-IFIX<DX*TRFACX^. 5)+KMINSX 

IYaIFIX(DY#TRFACY+. 5)-^KMINSY 

TREALX«=CX1 

TREALY-CYl 

CALL XYCNVTdX, lY) 

Et JD IF 

C ♦ 

C * PROCESS END POINT (GX2, CY2) 

C * 

BUFIN(l>eCXI 

BUFIN(2)='0Y1 

BUFIN(3>«*0X2 

BUFIN(4)«0Y2 


CniG.'NAL P,-G3 

VEC SUBROUTINE; Of POOR QUALITY 


VEC: VECTOR DRAW ROUTINE Cary Jonos/CSFC 

Oacenb»r 16, 1901 

PLOTIO requires two calls per vector, 'CALL MOVEA(GXl, CYl) ' 
and 'CALL DRAWA (CX2,GY2)'. This routine coflibines those 
operations into just one call with about a 20 percent reducti 
in run time. The calls in this subroutine are to PLOTIO. 
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CALL CLIPKDUFIN. BFOUT) 

IF (‘■'.CNFLC. NF. 1 > THEN 

modeSSde™^'"-*' 

CALL BUFFPK( 1,29) 

KMOVEF-1 

DX»CXS-TMINVX 

DYaCVS-TMINVY 

IX»IFIX{DX*TRFACX+. 5)-»KMlNSX 
IY=IFIX{DY*TnPACY-». 3)-»KMINFY 
CALL XYCNVT( IX. lY) 
KKMQDE-MODE 
END IF 

DX«CXE-TK1NVX 

DY-CVE-TMINVY 

IX=IFIX(DX*TRFACX+. 3)-»KI11NSX 
IY=»’FIX(DY*'*'RFACY-». 3)+KMINSY 
CALL XYCNVn IX. lY) 

TREALX--CXE 
TREALY^CYE 
END IF 

RETURN 

HrjD 







